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ABSTRACT 


A  variationally  consistent  tenth-order  displacement  theory  of  stretching  and  bending 
of  orthotropic  elastic  plates  is  proposed  which  lends  itself  perfectly  to  finite  element  for¬ 
mulations  based  upon  C°  and  C  '-continuous  displacement  approximations.  The  deforma¬ 
tions  due  to  all  strain  and  stress  components  are  accounted.  The  theory  is  derived  from 
three-dimensional  elasticity  via  the  principle  of  virtual  work  by  expanding  the  displace¬ 
ment  components  with  respect  to  the  thickness  coordinate  by  means  of  Legendre  polyno¬ 
mials,  where  the  transverse  displacement  is  of  a  special  parabolic  form  while  the  inplanc 
displacements  are  linear.  The  issues  of  thickness-expansion  related  inconsistencies  in  the 
transverse  shear  strains  and  the  transverse  normal  stress  are  resolved  in  a  rational  fash¬ 
ion.  The  resulting  parabolic  shear  strains  incorporate  Reissner’s  shear  correction  factor, 
while  the  transverse  normal  strain  varies  cubically  across  the  plate  thickness.  The  varia¬ 
tional  principle  yields  seven  equations  of  motion  and  e.xclusively  Poisson-type  edge  bound 
ary  conditions.  .A  qualitative  assessment  of  the  theory  is  carried  out  for  the  problem  of 
static  equilibrium  involving  an  infinite  plate  under  a  sinusoidal  normal  pressure.  Perti¬ 
nent  issues  on  the  particular  suitability  of  the  theory  for  the  development  of  efficient 
displacement  plate  elements  are  discussed. 
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INTRODUCTION 


In  many  design-critical  regions  in  plate  structures,  such  as  those  involving  highly  stressed 
areas  near  holes,  cutouts  and  impact,  bolted  and  adhesive  joints,  etc.,  deformations  due  to 
transverse  shear  and  transverse  normal  strains  and  stiesses  can  significantly  contribute  to  the 
three-dimensional  stress  state.  These  effects  are  especially  pronounced  in  composite  laminates 
that  are  inherently  compliant  in  the  traiuver.se  shear  and  transverse  normal  modes  of  deforma¬ 
tion.  To  adequately  model  such  critical  plate  regions  with  finite  elements,  three-dimensional 
elements  are  commonly  used,  usually  at  a  very  high  computational  cost.  Alternatively,  from 
the  standpoint  of  computational  efficiency,  and  in  order  to  avoid  or  minimize  costly  sub¬ 
structuring  with  three-dimensional  elements,  it  would  be  most  desirable  to  discretize  with  accu¬ 
rate  plate  elements  that  incorporate  all  strain  and  strc.s.s  effects. 

Many  effective  plate  finite  elements  developed  to  date  evolved  from  a  shear-deformable 
displacement  theory  originally  derived  for  homogeneous  isotropic  plates  by  Mindlin.’  The  the¬ 
ory  lends  itself  to  the  development  of  simple  and  efficient  elements  by  virtue  of  the  inherent 
C"-continuity  prescribed  for  the  approximation  of  the  plate  displacement  varhihlcs.  An  earlier 
plate  theory  of  Reissner""^  is  based  on  stress  approximations  including  the  transverse  normal 
stress  component  (which  is  neglected  in  the  Mindlin  theory)  and  is  somewhat  more  accurate 
as  far  as  the  stresses  are  concerned.  Both  of  the  Mindlin  and  Reissner  theories  rely  upon 
weighted-average  displacement  variables  which  constitute  linear  variations  across  the  thickness 
for  the  inplanc  displacements  and  a  uniform  tran.sverse  displacement.  From  a  perspective  of 
composite  laminate  analysis,  the  displacement-based  theories  arc  generally  preferred,  allowing 
for  discontinuous  normal  stresses  which  occur  naturally  in  composite  laminates  (e.g.,  refer  to 
Reference  5). 

In  recent  years,  several  higher-order  shear  deformable  displacement  plate  theories  have 
been  introduced  (e.g.,  see  References  6  and  7),  in  which  the  inplanc  displacements  vary  cubi- 
cally  across  the  plate  thickness  while  the  transverse  displacement  is  uniform.  While  these  the¬ 
ories  appear  to  have  a  wider  applicability  range  than  their  predecessors,'’^  they  possess  a 
serious  drawback  from  the  finite  clement  standpoint,  the  requirement  of  C'  continuity  which 
is  known  to  inhibit  development  of  simple  and  efficient  elements  (e.g.,  sec  discussions  in 
Reference  8). 


Higher-order  plate  theories, which  in  addition  to  transverse  shear  also  include  the 
effect  of  transverse  normal  straining,  arc  intended  to  extend  the  range  of  applicability  to 
design-critical  situations  such  as  those  mentioned  earlier.  The  interested  reader  may  find 
insightful  reviews  of  many  notable  higher-order  theories  in  the  papers  by  Lo,  Christensen,  and 
Wu  "^  and  Reissner.'"  Regrettably,  these  theories  are  rather  unattractive  from  the  finite  ele¬ 
ment  viewpoint  for  the  following  reasons:  (i)  an  excessively  large  number  of  plate  displace¬ 
ment  variables  (e.g.,  Lo  ct  al.’s  theory'^  possesses  11  kinematic  variables  for  the  problem  of 
bending  and  stretching,  whereas  the  Reissner  and  Mindlin  theories  involve  only  five  such  vari¬ 
ables),  (ii)  natural  edge-boundary  conditions  that  include  nonclassical  quantities  such  as  higher- 
order  force  and  moment  resultants'^’*^  (in  addition  to  the  physically  meaningful  Poisson 
conditions  including  three  prescribed  edge  forces,  bending,  and  twisting  moments),  (iii)  a  lack 
of  variational  basis, and  (iv)  a  departure  from  the  kinematically  admissible  evaluation  of 
stress;'^  i.e.,  evaluation  of  stre.ss  co.mponents  from  equilibrium  equations  rather  than  (in  a  kine¬ 
matically  admissible  manner)  from  constitutive  relations. 
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The  present  development  of  a  higher-order  plate  theory,  incorporating  all  stress  and  strain 
effects,  is  motivated  by  finite  element  application,  where  a  variational  formulation  is  desired 
requiring,  at  most,  C’-continuous  displacement  fields  with  the  least  number  of  variables,  and 
relying  exclusively  upon  Poisson-type  boundary  conditions. 

In  the  Kinematic  Assumptions  Section,  we  expand  three  Cartesian  displacements  in  terms 
of  a  thickness  coordinate  using  Legendre  polynomials  to  simplify  integration  across  the  thick¬ 
ness.  To  obtain  the  simplest  higher-order  theory,  we  adopi  the  same  expansion  order  as  pro¬ 
posed  by  Hildebrand,  Reissner,  and  Thomas;^  linear  and  parabolic  distributions  for  the  inplane 
and  transverse  displacements,  respectively.  In  contrast  to  Hildebrand  et  al.’s  transverse  dis¬ 
placement  expansion  in  which  the  leading  term  is  the  midplane  displacement,  the  leading  term 
in  the  present  expansion  is  Reissner’s'^  weighted-average  displacement  Here,  as  in  Ref¬ 
erence  9  the  parabolic  transverse  displacement  also  involves  two  higher-order  dispiacemenl 
components.  The  particular  advantage  of  the  present  expansion  strategy  is  manifested  in 
finite  clement  approximations,  where  all  kinematic  variables  need  not  exceed  C°  continuity. 

In  the  Concluding  Summary  and  Finite  Element  Suitability  Section,  we  shall  remark  further 
on  some  important  consequences  of  the  pre.sent  continuity  requirements  on  finite  element 
development. 

In  the  Transverse  Shear  Field  Consistency  and  Transverse  N'ormai  Field  Consistency  Sec¬ 
tions,  we  delineate  certain  types  of  inconsistencies  that  arise  in  the  thickness  expansion  for 
the  transverse  shear  strains  and  transverse  normal  stress,  and  propose  a  novel  approach  which 
eradicates  these  inconsistencies  in  a  rational  manner.  In  the  Governing  Equations  Section,  we 
employ  the  principle  of  virtual  work,  which  yields  seven  partial  differential  equations  of 
motion  and  exclusively  Poisson  edge  boundary  conditions.  In  the  special  case  of  static  equilib¬ 
rium,  the  theory  decouples  into  one  of  fourth-order  inplane  stretching,  sixth-order  transverse 
bending,  and  zero-order  transverse  stretching.  Whereas  the  first  two  sets  of  equations  corre¬ 
spond  to  Reissner’s  theory, the  last  set  includes  two  higher-order  transverse  displacements 
of  zero-order  (i.e.,  they  possess  no  spatial  derivatives)  and  first-order  derivatives  of  the 
weighted-average  variables.  The  solution  to  a  plate  boundary  value  problem  simply  involves 
solving  the  fourth-order  inplane  stretching  and  sixth-order  transverse  bending  equations  subject 
to  Poisson  boundary  conditions  (as  in  Reissner’s  theory),  and  then  a  straightforward  evaluation 
of  the  two  higher-order  variables  from  the  zero-order  transverse  stretching  equations. 

One  remarkable  aspect  of  the  theory  is  the  particularly  accurate  representation  of  all 
transverse  stresses.  Thus,  the  kinematically  admissible  transverse  shear  and  normal  .stresses 
are  shown  to  satisfy  exactly  the  transverse  equilibrium  equation  of  elasticity,  i,e.,  thev  arc  also 
statically  admissible  in  the  transverse  equilibrium  ,scn.sc.  The  transverse  shear  stresses  vary  par- 
abolically  across  the  thickness  while  the  transverse  normal  stress  has  a  cubic  distribution.  .-\11 
transverse  stresses  satisfy  the  prescribed  tractions  at  the  top  and  bottom  plate  faces. 

In  the  Assessment  of  Theory  Section,  we  examine  the  range  of  applicability  of  the  theory 
via  an  analytic  solution  for  an  infinite  isotropic  plate  under  quasi-static  sinusoidal  normal  pres¬ 
sure  and  compare  results  with  the  exact  elasticity  solution  and  several  other  plate  theories. 

Finally,  in  the  Qtncluding  Summary  and  Finite  Element  Suitability  Section,  we  discuss  the 
issues  concerning  utilization  of  this  variational  displacement  theory  within  the  finite  element 
framework,  and  suggest  ways  of  formulating  displacement  plate  elements  whose  computaticmal 
efficiency  would  be  com.parable  to  Reissner-Mindlin  type  elements.'*' 
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KINEMATIC  ASSUMPTIONS 


Consider  an  orthotropic,  homogeneous  elastic  plate  located  in  the  three-dimensional  Carte¬ 
sian  framework  (x,  y,  z),  where  the  x-y  plane  passes  thiough  the  middle  of  the  plate  thickness 
I  =  0.  where  |  =  z/h  e[-l,  1]  is  a  dimensionless  thickness  coordinate,  and  2h  is  the  plate 
thickness.  The  tun  (S^)  and  bottom  (S’)  plate  surfaces  are  taken  free  of  shear  tractions  and 
loaded  by  normal  pressures 


(i=x,y). 

=  q^(x,y,t) 

on 

(1) 

(i=x,y) , 

a  (5=-l)=  q  (x,y,t) 

on  S 

r-) 

where  t  denotes  time.  On  a  part,  S^j,  of  the  cylindrical  boundary  surface  S  normal  to  the 
midplane  Sjn,  a  traction  vector  is  prescribed: 


T  ,  1  }  on  £ 
y  2  o 


U’e  as.sume  the  plate  principal  material  directions  to  be  coincidents  with  the  Cartesian 
directions,  and  small  displacements. 


(3) 


The  simplest  physically  meaningful  approximations  of  displacement  components  allowing 
for  the  transverse  normal  straining  involve  linear  through-thicknes.'  expansions  of  the  inplane 
displacements  (as  in  classical  and  refined  theories),  and  a  parabolic  approximation  of  the  trans¬ 
verse  displacement.  This  level  of  approximation  was  first  discussed  by  Hildebrand  et  al.^  For 
convenience,  we  expand  the  Cartesian  displacement  components  by  means  of  Legendre  polyno¬ 
mials  with  Uz  having  a  special  form; 


u  (x.y.z.t)  =  u(x,y,t)  +  hPi(C)6,^(x,y,t) , 

X  y 

u  (x,y,z,t)  =  >/(x,y,t)  +  hPi(5)0„(x,y,  t), 

y  ^ 

u^(x,y,z,t)  =  w(x,y,t)  +  Pi(c)wi(x,y,t)  +  [ Po (C )/5+P2 ( C) jw- (x, y , t ) , 

wncre 

p^'-;)  =  — —  ?  (i)=i,  (5=  z/h  c[-i,i]) 

d""  n!  ^ 

(i.e.,  P3=  1,  P:=  P;=  (3^-U/2,  P3=  4(5^'-3)/2  ) 

are  Legendre  polynomiaN  which,  due  to  their  orthogonality  property, 
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P  P 
m  n 


-1 


0  if  m?=n 

2 / ( 2m+ I )  if  m=n 


(^) 


prove  convenient  whenever  thickness  integrations  are  performed. 

The  expansion  coefficients  u,  v,  6^,  0^,  and  w  in  Equation  4  are  weighted  average 
kinematic  variables  defined  here  in  accordance  with  Reissner.'^ 
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X  2h^  J 
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u  z  dz  ,  0 
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u  z  dz 
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-h 


3 

2h^ 


h 

u^(P(3-  P2)  dz,  {  Pq 
-h 


P2  =  3(1  -  5')/2  }, 


(2) 


(8) 


(9) 


where  u  and  v  denote  the  midplane  displacements  along  the  x  and  y  directions,  respectively 
and  and  dy  denote  the  normal  rotations  about  the  x  and  y  axes,  respectively.  In  Equa¬ 
tion  4,  w,  and  w;  are  the  higher-order  components  of  a  transverse  displacement,  u^,  which 
varies  parabolically  across  the  thickness.  Note  that  w  is  the  weighted  average  transverse  dis¬ 
placement,  in  contrast  to  ether  higher-order  theories^''^  employing  the  midplane  displacement 
variable.  As  will  be  seen,  the  special  expansion  for  u^  leads  to  exclusively  Poisson-type  edge 
boundary  conditions  which  are  naturally  obtained  from  the  virtual  work  principle. 


The  plate  strain-displacement  relations  are  obtained  in  the  usual  manner  by  introducing 
Equation  4  into  the  strain-displacement  relations  of  three-dimensional  elasticity; 


e .  =  u. , . , 
1  11 


''lj“  ‘"i’j  “j’i'’ 


resulting  in 


e.  =  u,  +  ,  ,  e  ■=  v,  +  hC9  , 

X  ’x  ^  y’x’  y  y  x’y 


f  =  u,  +  V ,  +  h5(8  ,  +  9  ,  ) 

XV  y  X  ^  V  V  XX 


(10) 


(11a) 
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3 

0 
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‘  X 

1 

3 

3 

(lib 

'  To 

-t-  + 

2 

^  '  V 

Y 

8  +  w 

.  >  Y 

=  9  +  w. 

XZq 

y 

X  'yzo 

X  y 

=  (w;  -t-  2E,w2)/h  (11c) 

where  a  comma  (,)  denotes  partial  differentiation  with  respect  to  the  spatial  coordinates.  To 
facilitate  our  subsequent  discussion,  we  expressed  the  strains  in  Equations  11  directly  in  terms 
of  the  thickness  coordinate  q. 
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The  complete  Hooke's  law  for  a  homogeneous  orthotropic  plate  may  be  expressed  in 
matrix  form  as; 

O  =  C  E, 
or 
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(Symm. ) 
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TRANSVERSE  SHEAR  FIELD  CONSISTENCY 

From  Equation  11b,  it  becomes  evident  that  certain  ^-distribution  inconsistencies  are  pres¬ 
ent  in  the  expansions  of  '/xz  and  yjj_.  Here,  both  linear  and  quadratic  terms  in  |  involve  only 
wi  and  W2,  which  are  contributed  from  u^.  This  is  obviously  the  result  of  Uj,^  (i  =  x,  y) 
being  uniform  across  the  thickness,  whereas  vary  parabolically.  It  is  clear  that  to  possess 
a  consistent  kinematic  form,  each  ^-polynomial  coefficient  must  include  contributions  from 
both  the  inplane  and  transverse  displacements,  as  implied  in  Equation  10. 

A  simple  way  of  assessing  the  effect  of  the  aforementioned  inconsistency  is  to  examine 
the  nature  of  the  transverse  shear  strains  for  a  plate  under  the  Kirchhoff  constraint,  i.e., 

(i  “  y)  the  plate  thickness  diminishing  to  zero.  From  each  shear  strain,  there 
results  three  constraint  equations; 

J  -  0 

t  i i 

j  -  0  ■ 

Assessment  of  Theory  Section); 


{  V 

'XZo 

To 

V  ^ 

^  ’x’ 

'Vo 

3 

10 

-  ’  y’ 

"-y 

It  can  be  shown  that  (e.g.,  see  results  in  the 


=  oClh/D" 

^  iz^ 


-^  =  0f2h/L)S  (i=x,7) 
^  iz. 


(13a) 


a 


where  L  denotes  some  characteristic  plate  span.  .Accounting  for  Equation  13a  in  Equa¬ 
tion  13,  it  becomes  apparent  that,  at  least  in  the  Kirchhoff  regime,  the  resulting  shear  strains 
are  uniform.  css  the  plate  thickness; 

-*0,  ( i=x,y)  .  (14) 

The  direct  im.plication  of  Equation  14  is  that  explicit  conditions  of  zero  shear  tractions  on  the 

plate  faces  cannot  be  satisfied  and,  moreover,  an  ad  hoc  specification  of  a  shear  correction 

factor  is  necessarv  to  achieve  agreem.ent  with  the  Kirchhoff  theorv  for  the  class  of  thin  plate 
,,  !!,i:  -  s  -  I 

problems. 

To  eradicate  this  type  of  kinematic  inconsistency  while  retaining  a  parabolic  distribution 
for  u^,  the  gradients  u.^  (i  =  x.  y)  should  naturally  be  represented  by  a  complete  quadratic 
polynomial  in  g,  thus  implying  cubic  distributions  for  u,  and  Uy.  Approximations  of  this  type 
were  employed  by  Lo  et  al.^  (who  extended  Reissner’s  work'"*  in  which  inplane  modes  of 
deformation  were  neglected)  resulting  in  a  complex  twenty-second-ordcr  theory  of  limited  prac¬ 
tical  usefulness. 


In  the  present  treatment,  in  order  to  resolve  this  difficulty  while  preserving  the  simpliciiN 
of  Equation  -i,  we  propose  to  replace  the  uniformly  distributed  gradients  u,j  (i  =  x.  y; 
(directly  computed  from  Equations  4  and  10)  with  the  field-consistent  gradients  u*,^  v'.hicn 
vary  parabolically  across  the  thickness: 


2 

u.,^  =  ^  a.j^(x,y,t)  (i=x,y). 

k=0 


(15) 


Henceforth,  the  superior  asterisk  (*)  will  denote  field-consistent  gradient  and  resulting  stress 
quantities.  The  relations  between  a,|;  and  plate  displacement  variables  are  determined  via  two 
types  of  physically  desirable  constraint  conditions: 

The  field-consistent  transverse  shear  stresses  are  required  to  satisfy  tract. on-free  boundary  con¬ 
ditions  at  the  t(5p  and  bottom  plate  faces; 


=  0  ( i=x, y) .  ( lu) 

(5=1,  -•-) 


Taking  into  account  the  constitutive  relations.  Equation  12.  the  above  stress  constraints  give 
rise  to  the  vanishing  conditions  upon  the  shear  strair*': 


Y  •  =  { u .  ,  -I-  u  ,  .  } 

iz  1  z  z  1 


=  0  (  i=x,y) , 

(5=1,  -1) 


t  Ida  i 


Furthcrmoirc.  the  field-consistent  shear  strains  arc  required  to  be  equivalent  in  the  mean  to 
their  field-inconsistent  counterparts  from  Equation  11b: 
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min  U  .  -  Y  .  )■  dz,  ( i=x,y). 

J  iz  iz 


Equations  I6a  and  17  are  solved  I'or  a,k  (i  =  x,  y;  k 
shear  strains: 


=  0,  I,  2),  giving  rise  to  the  consistent 


in  which  k.-  =  5,t)  emerged  from  this  derivation  and  is  coincident  with  Rcissner’s  shear 
correction  factor. 

TRANSVERSE  NORMAL  FIELD  CONSISTENCY 

.A  stress  type  of  field  inconsistency  can  be  aetccted  in  the  transverse  equilibrium  equation 
of  the  three-dimensional  elasticity; 

T  -*■  T  +  a  =0  (ignoring  the  body  force)  (19 

xz,x  yz.y  z.z 

Here,  owing  to  the  parabolic  g-distribution  of  r,^  (i  =  x,  y).  must  be  cubic  in  however, 

assumptions,  Equation  4,  can  only  yield  a  which  is  linear.  To  overcome  this  discrepancy,  \\'c 
replace  u^^  (computed  directly  from  Equations  4  and  10)  with  the  field-consistent  gradient 
u*2,2,  having  a  cubic  ^-distribution  across  the  thickness: 

j 

q  -  -  I  t,yy.,y,t)  P^(0  CO 


vncre  b,^  .ire  determined  irom  the  :..:dc/Aing  constraint  equations: 


The  homoecr'.eous  bound.irv  conditions  ’..in  the  trtrough-thickness  gradient  oi  transverse  norm.'.i 
stress  component: 

»'  =0,.  : 

These  conditions  result  from  Equation  19  when  taking  the  top  and  bottom  plate  faces  to  be 
free  of  shear  tractitins. 

The  field-consistent  transserse  normal  strain  is  required  to  be  equivalent  in  the  mean  to  its 
fie;d-inc:ins:,stcnt  counterpart  ot  Equation  lie; 


min  .  <  -s 


Iz  . 


F^'orn  Equations  21  and  22,  there  results  tour  equations  in  terms  of  b^;  yielding  a  cubic 
transverse  normal  strain  of  the  form; 


with 
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+  k' 
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L (bP^-Pj )w; 


h- 
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(UP3+  P.)(C 
3  3 


3  L 


d  + 
y.x 


(23) 


k‘  =  42/85 
2 

In  the  remaining  part  of  the  formulation  we  shall  exclusively  deal  with  the  field-consistent 
strains  and  stres.ses  and,  hence,  we  shall  omit  the  superior  asterisk  which  was  previously  used 
to  distinguish  these  quantities. 


GOVERNING  EQUATIONS 

We  derive  the  equations  of  motion  together  with  the  natural  boundary  conditions  via  the 
principle  of  virtual  work  incorporating  the  tield-consistent  strains  and  stresses: 


t  ; 

r  '  ' 

i  <  i  !  (o  0£  +0  -Sc  ~  <5  5c  T  5y  +  t  OV  +  t  5-^  )  dxdvc 

j  i,  J  J  J  y  y  ^  ^  xy  xy  X2  X2  y2  vz 

t-  V 


4  0  (u  +  u  +  u  )  dxdydz 
X  y  z 


(24) 


-  f  q^ou  dxdy  - 

Jj  z 


q  5u^dxdy 


(T  6u  +  T  5u  +  T  5u  )  dsdz  ^  'dt=  0 
X  X  y  y  z  z 


where  s  is  a  coordinate  tangent  to  the  boundary  of  S^j,  p  is  the  mass  density,  and  the 

superior  dot  (■)  denotes  differentiation  with  respect  to  time. 

Consistent  with  the  parabolic  distribution  of  the  shear  strains.  Equation  13,  the  prcscrilted 
edge  shear  traction  is  assumed  to  vary  parabolicaily  as  well,  i.e.,  T__  =  (P)  -  P;)  T-(s).  Intro¬ 
ducing  Equations  l!a.  1?  and  23  into  Equation  24.  and  integrating  with  respect  t.a  z,  we  find 

that  certain  intearals  vanish  identicailv 


!  T  ?  ■  dz  =  '1  T  (Pi/5  -  ?-}  Qz  =  a  (,P.+14p,j  dz  =  0, 
i  2  •  j  z  ■  ‘  J  z  ‘  ^ 

-h  -h  -h 


'  2o  a ! 


yielding  a  two-dimensional  virtual  work  statement  of  the  form; 

t:  r 

iff  ^ 
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Performing  appropriate  variations  and  integration  by  parts  yields; 
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where  C,  and  C,j  are  parts  of  the  boundary  surrounding  S^i  where  the  tractions  and  displace¬ 
ments  are  prescribed,  respectively;  and  where  the  inertial  and  stress  resultants  arc  defined  as; 
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(N  ,  N  )  =  (T  ,  T  )  dz,  (M  ,  M  )  = 
xn  yn  J  x’  y  xn’  yn 
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zn 


T  dz, 
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N  =  N  i  +  N  ra,  N  =  N  £  +  N  m, 

xn  X  xy  yn  xy  y 

M  =  Ml+M  m,  M  =M  l+Mm,  Q  =  Q£  +  Qm, 

xn  x  xy  yn  xy  y  zn  x  y 

2,  =  cos(x,n),  m  =  cos(y,n), 

where  n  is  an  outward  normal  to  the  boundary’  cl'  Sqj. 

Thus,  the  principle  provides  seven  equations  of  motion; 
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( Sw; ) : 


-  M  /  h  +  q  •  /  k  ■ 
z 


=  :  (w  +  2W2) 


•Jr’) 


five  stress  boundarv'  conditions  on  C,/ 
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xn  xn 


N  =  N  , 
■m  vn 


M  =  M  , 
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(29a) 


and  five  displacement  boundary  conditions  on  C^; 


u=u,  v=v,  0=9,  9=0,  w=w. 

y  y  X  X 


(29b) 


Note  that  Equation  29  represents  Poisson-type  boundary  conditions. 

Invoking  Hooke’s  law  for  an  orthotropic  material.  Equation  12.  yields  the  plate  stress- 
displacement  relations  of  the  form: 
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A.  .  =  (2h)  C.  . ,  (i,j  =  1.2,3), 

i  j 
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ij 


C.  .  - 
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( i, j=l,2). 
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Special  Case:  Plate  in  Equilibrium 


It  is  of  significance  to  consider  the  case  of  a  quasi-static  loading,  leading  to  the  vanishing 
of  the  inertial  terms,  in  which  case  Equations  28a  through  28g  become  equilibrium  equations. 
Substituting  Equations  30i  and  30j  into  Equations  28f  and  28g,  respectively,  results  in  the 
solution  for  wi  and  W2: 


w , 
h" 


'Z2  - 
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‘3  3 


qjh  - 


q,  >2  ■  D139  1  '  . 

i3y>x  ■‘^xy 


(32a) 
(32  b  I 


Introducing  Equation  32  into  Equation  30  yields  the  expressions  for  Nx,  Ny,  Mx,  and  Myi 

0 1 382ft 
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where 
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When  Equations  30c,  30f  through  30h,  and  33  are  introduced  into  the  equilibrium  Equations 
28a  through  28e,  we  obtain  two  decoupled  sets  of  second-order  partial  differential  equations 
of  equilibrium  in  terms  of  five  weighted-average  displacements: 

Fourth-order  stretching 


A,,u,  -t-  A12V,  +  Afifi(u,  +  V,  )  + 

«xx  xy  yy  ’xy 


+  A22V,  +  AftfiCu,  +  V,  )  + 

w  xy  yy  xy  xx 
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^33  ^ 
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(35a) 
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Sixth-order  bending 


,  -t-  Di.9  ,  +  DggO  ,  +  9  ,  )  -  GssCw.  8  )  =  0 
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’xx  y’x  yy  x’y 


(35  b) 


with  these  equilibrium  equations  being  subject  to  the  natural  boundary  conditions.  Equation  29. 

*■  .  ^  ■'-4 

It  is  quite  remarkable  that  comparing  these  results  to  Reissner’s  theory”  (derived  lor  the 
case  of  transverse  bending  of  isotropic  plates,  i.e.,  excluding  the  inplane  deformations  u  and 
v),  we  find  that  the  equilibrium  Equation  35b,  the  corresponding  boundary  conditions.  Equa¬ 
tion  29,  and  the  stress  resultants.  Equations  33b,  and  30f  through  30h  match  those  ot  the 
Reissner  theory. 


Computing  stresses  from  the  constitutive  relations.  Equation  12,  we  find  a  cubic  trans¬ 
verse  normal  stress 


a 


z 


1 

4 


qiC(3 


C')  +  2q2 


(36) 


satisfying  the  prescribed  normal  pressure  boundary  conditions  of  Equations  1  and  2.  From 
Equation  12,  we  also  determine  the  parabolic  shear  stresses  of  the  form 


3k^ 

C4,(l 

-  ^M(w,  +•  0  ) 

= 

Q  (1 
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.X  y 

4’n 

X 

- 

c.,( : 
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_  3_ 

■  yz  2 

V  X 
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These  kinematically  admissible  t'-ansverse  stresses  (i.e.,  obtained  from  constitutive  equations) 
can  be  shown  to  satisfy  the  transverse  equilibrium  equation  of  elasticity.  Equation  19,  and. 
hence,  they  are  also  statically  admissible.  The  transverse  stresses  can  also  be  seen  coincident 
with  those  of  Reissner,”"'  who  assumed  the  parabolic  transverse  shear  stresses  at  the  outset, 
and  then  computed  the  transverse  normal  stress  from  Equation  19.  As  far  as  the  inplane 
stresses,  the  present  identical  to  Reissner’s,  however,  and  Oy  are  cubic  across  the 

thickness  due  to  Equation  23,  whereas  in  Reissner’s  theory  they  are  linear. 

Note  that  although  in  statics  the  equations  of  equilibrium  are  conveniently  decoupled,  in 
dynamics,  all  seven  equations  of  motion.  Equation  28,  are  fully  coupled.  In  the  latter  case, 
the  solution  procedure  involves  substitution  of  the  stress  resultants  Equation  30  into  Equa¬ 
tion  28,  and  simultaneously  solving  the  resulting  seven  differential  equations  for  the  kinematic 
unknowns. 
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ASSESSMENT  OF  THEORY 

To  quaiitati -elv  evaluate  the  present  plate  theory,  an  analytic  solution  is  obtained  tor  the 
problem  of  an  infinite  plate  in  equilibrium  loaded  at  the  top  surface  by  a  sinusoidal  pressure 
of  the  form: 

q"*"  -  qo  sin(TTx/L)  on  S  (q  =  0  on  S  )  (38) 


where  L  is  a  half  wavelength  of  the  loading  pattern.  This  problem,  having  an  exact  elasticity 
solution,*^  has  been  widely  employed  as  a  benchmark  in  assessing  plate  bending  theories. 


Introducing  Equation  38  into  Equation  35,  and  seeking  the  unknown  kinematic  variables 
in  the  form  of  sin^ntx/L)  and  cos(--rx/L),  a  system  of  ordinary  differential  equations  is  obtained 
and  readily  solved  for  the  unknowns.  The  solutions  are  then  simplified  for  the  isotropic  case 
by  replacing  QjS  from  Equation  12  with  the  appropriate  expressions  in  terms  of  the  Lame' 
constants,  /  and  u,  which  result  in  the  form: 


v(x)  =  0  (x)  =  0 

X 


u(x) 


SiTuC  A+u) 


COs(  'TX/  L  j 


e^(x) 


5qoL^ 


A[1  +  2u/A  -  (2/5)(iTh/L)^] 
lj(A+y) 


cosCirx/L) 


(39) 


..  ,'..s  _  3qoL‘‘  A[1  +  2u/A  +  ( 2/ 5 ) (rrh/L) ^( 3  +  4y/A)] 

-  8^  - ilTA+^T) - ^  sinCnx/L)  . 

Relations,  Equation  32,  are  then  invoked  to  determine  the  remaining  kinematic  variables: 


sin(Trx/L) 


8u(X+u) 

.  +  2u/\  -  (2/5)(-Th/L)Ml  *  85m/21A)]  ,  ,  ,,  , 

uJf.)  =  -  - ^ -vr'-o  \ -  sin(TTx/L) 

‘  86-^  n  yt  A+M  j  I  at-2u) 


(40) 


The  midplane  transverse  displacement  is  obtained  from  Equation  4: 


u  (x,5=0)  =  w(x)  -  0.3  w-(x) 
z 


(  4  I  ; 


Hooke’s  law,  Equation  12,  then  provides  a  means  for  obtaining  stresses: 

sin(TTx/L) 


6(L/.h)«+  ,  J  -  5  =  ) 

A  +  2y  5 


o  (x,^)  =  ^  (2-0(C+f)^  sin(iTx/L) 

z  A 


(42) 


:  (x,^)  =  cos(-irx/L) 

xz  Atth 
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Several  observations  can  be  made  regarding  these  results: 

(a)  .As  expected,  all  weighted-average  kinematic  variables  are  identical  to  those  of  the 
Reissner  theorv'.  One  can,  therefore,  regard  the  expansion  for  as  hierarchical,  in  which 
Reissnei’s  displacement  w  is  the  leading  term,  and  Wj  and  wt  are  the  higher-order  terms  allow¬ 
ing  a  parabolic  variation  of  across  the  plate  thickness. 

(b)  In  the  axial  stress  component  a.^,  the  first  term  (linear  in  .;)  is  the  Reissner  (or  classi¬ 
cal)  stress,  while  the  second  term  is  a  higher-order  contribution  varying  as  Pi(g).  .Vote  that 
the  leading  stress,  of  the  order  (L/h)',  strongly  dominates  with  the  decreasing  plate  thickness. 

(c)  .As  established  earlier,  and  are  the  same  as  in  the  Reissner  theory;  they  satisfy 
both  constitutive  equations  and  the  transverse  equilibrium  equation  of  elasticity. 

In  Figure  1,  a  normalized  value  of  the  maximum  midplane  transverse  displacement  is  plot¬ 
ted  versus  the  2h,l  ratio,  where  the  results  of  five  bending  theories  arc  compared  with  the 
exact  solution.  It  is  seen  that  the  present  and  Essenburg’s'^  solutions  agree  quite  closelv, 
and  they  begin  to  deviate  from  the  exact  solution  at  2h,^L  =  1.  This  is  anticipated  since  the 
two  theories  use  the  same  order  of  displacement  approximations,  the  major  difference  being 
that  in  Essenberg's  beam  theory,  both  stresses  and  displacements  are  assumed  indepencentiy. 
The  eleven-variable  theory  of  Lo  et  al.'^  also  appears  to  be  accurate  up  to  the  same  ratio 
value.  These  esults  indicate  that  when  a  characteristic  length  of  loading  (or  plate  span)  is 
of  the  order  of  the  plate  thickness,  application  of  a  higher-order  theory  is  particularly  war¬ 
ranted.  Similar  observations  have  been  made  with  regard  to  contact  problems.’"^  as  w'cll  as 
those  involving  disturbances  of  geometry  such  as  holes  and  cutouts,  having  the  characteristic 
length  of  the  order  of  the  plate  thickness.  Note  that  even  though  Reissner’s  weighted- 
average  displacement,  w,  is  also  shown  in  the  comparison  (the  only  transverse  displacement 
variable  in  Reissner  theory),  it  does  not  represent  the  midplane  displacement.  In  contrast,  w 
is  the  leading  term  in  the  expansion  of  the  present  theory,  where  the  complete  through-the- 
thickness  distribution  is  expressed  by  u^  in  Equation  4. 

-  Exact 


C«l/  l_ 


Figure  1.  Maximum  midplane  transverse  displacement 
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In  Figures  2  through  4,  the  present  results  for  a^,  and  are  presented  for  the  case 
of  2h.L  =  0.5,  and  compared  with  the  corresponding  exact  stresses.  It  is  seen  that  the  integ¬ 
rity  of  all  stresses  is  excellent,  and  particularly  the  transverse  shear  and  normal  stresses  for 
which  other  displacement  theories  produce  rather  inadequate  stresses  violating  stress  boundary 
conditions  on  the  top  and  bottom  plate  faces  (e.g..  References  12  and  15). 


Figure  2.  Distribution  of  maximum  ox 
stress  across  thickness, 

CONCLUDING  SUMMARY  AND  FINITE  ELEMENT  SUITABILITY 

.‘X  tenth-(?rder  displacement  theory  for  orthotropic,  homogeneous  elastic  plates  which 
includes  the  inplane,  transverse  shear,  and  transverse  normal  deformations  was  developed. 

The  theory  embodies  two  more  kinematic  variables  than  the  Reissner  and  .Mindlin  shear- 
deformable  theories.  A  novel  means  for  removing  field  inconsistency  in  the  transverse  shear 
.strains  was  implemented  and  led  to  parabolic  strain  distributions,  satisfying  traction-free  bound¬ 
ary'  conditions.  In  a  similar  fashion,  the  transverse  normal  strain  and  stress  approximations 
were  raised  to  cubic  thickness  variations.  The  virtual  work  principle  yielded  seven  equations 
of  motion  together  with  exclusively  Poisson-type  boundary  conditions.  In  the  case  of  static 
equilibrium,  the  theory  revealed  a  close  relationship  to  Reissner’s  theory.  An  analytic  solu¬ 
tion  to  an  infinite  isotropic  plate  under  a  sinusoidal  load  revealed  the  high  order  of  accuracy 
attainable  with  this  theory  for  all  kinematic  and  stress  quantities  and,  particularly,  the  trans¬ 
verse  stresses  for  which  previous  displacement  theories  yielded  inaccurate  results. 
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Finally,  the  proposed  displacement  theory  appears  to  be  ideally  suited  for  finite  element 
approximation.  The  five  weighted-average  kinematic  variables  require  only  C“-continuity 
across  element  boundaries,  having  the  highest  spatial  derivative  in  Equation  25b  of  order  one. 
The  remaining  two  displacements,  w^  and  Wn,  do  not  possess  derivitives  in  the  energy  expres¬ 
sions,  thus  permitting  interelement  discontinuous  shape  functions  (i.e.,  C’  continuous).  The 
latter  aspect  allows  the  static  condensation  of  degrees  of  freedom  associated  with  wj  and  W2 
at  the  element  level.  The  implication  is  that  computationally  efficient  elements  can  be  devel¬ 
oped  by  directly  adopting  successful  approaches  of  shear-deformable  elements. Efforts 
are  currently  underway  to  develop  beam  and  plate  elements  based  on  this  higher-order  dis¬ 
placement  theory. 
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)  displacement  plate  elements  are  discussed. 


